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The PAMELA positron excess is analyzed within the framework of nonuniversal SUGRA models 
with an extended C/(l) n gauge symmetry in the hidden sector leading to neutralino dark matter with 
either a mixed Higgsino-wino LSP or an essentially pure wino dominated LSP. The Higgsino-wino 
LSP can produce the observed PAMELA positron excess and satisfy relic density constraints in the 
extended class of models due to a near degeneracy of the mass spectrum of the extended neutralino 
sector with the LSP mass. The simultaneous satisfaction of the WMAP relic density data and 
the PAMELA data is accomplished through a co-annihilation mechanism (Bc — mechanism), 
and leads to predictions of a neutralino and a chargino in the mass range (180-200) GeV as well 
as low lying sparticles accessible at colliders. We show that the models are consistent with the 
anti-proton constraints from PAMELA as well as the photon flux data from EGRET and FERMI- 
LAT. Predictions for the scalar neutralino proton cross section relevant for the direct detection of 
dark matter are also discussed and signatures at the LHC for these PAMELA inspired models are 
analyzed. It is shown that the mixed Higgsino-wino LSP model will be discoverable with as little as 
1 fb -1 of data and is thus a prime candidate for discovery in the low luminosity runs at the LHC. 
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I. INTRODUCTION 

The recent results from the PAMELA experiment [l[ 
which indicate a large excess in positron production in 
the galactic halo, along with a lack of any significant ex- 
cess in the anti-proton flux in the same experiment, have 
resulted in a large effort to understand the data. Further, 
the results from the FERM1-LAT experiment relating 
to the electron and photon fluxes emanating from the 
galaxy have also recently been reported and several in- 
terpretations have been put forth. These include particle 
physics models d i, H H, 0, H !, M, ED, as well 
as models where the source of the PAMELA excess is 
of astrophysical origin Here, we give an analysis of 
the positron excess in a supergravity (SUGRA) frame- 
work with the neutralino as the lightest R-parity odd 
supcrsymmetric particle (LSP), where we go beyond the 
parameter space of the minimal supergravity grand uni- 



fication model [14j . As pointed out in Ref. [L5|, annihi- 
lations of neutralinos into W + W~ and their subsequent 
decay can provide a positron excess. One outstanding is- 
sue in supersymmetric interpretations of the positron ex- 
cess has been that models with large neutralino annihila- 
tion rates into W + W~ states tend to imply a thermally- 
produced cold dark matter relic density which is as much 
as 100 times smaller than the amount indicated by the 
recent WMAP data [ijj]. This problem has been dis- 



cussed in the context of a non-thermal mechanism which 
is capable of generating the right amount of the relic 
density @, S] ■ 

Such non-thermal mechanisms could take 
place via moduli decay which can arise in various softly 
broken supersymmetric theories 17 1. 



In this paper we consider an alternative mechanism 
that can generate a relic density for the LSP neutralino 
which utilizes the idea of coannihilation of the LSP with a 
cluster of states in the hidden sector which are degenerate 
in mass with the LSP. Phenomena of this type can arise in 
a broad class of models including extended SUGRA mod- 
els with a Fayct-Iliopoulos term, models with kinetic mix- 
ings, models with the Stueckelberg mechanism, and in 
string and D-brane models. In many such models connec- 
tor fields exist which have non-trivial gauge transforma- 
tions under the hidden as well as under the visible sector 
gauge groups and allow for mixings of the LSP with the 
hidden sector gauginos and chiral fcrmions. The mixings 
between the visible and the hidden sector are typically 
constrained by the precision electroweak data. There are 
many examples of models of this type that can be con- 
structed, but here we will consider one specific concrete 
manifestation as a representative of this class of mod- 
els. More specifically, we will assume that there is a set 
of connector fields (axions) that transform non-trivially 
under the hidden sector U(\) x gauge group (i.e., the 
product U(\)x x U(l)x' • ■■) as well as under the hyper- 
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charge gauge group U{\)y- In addition to the above, in 
the visible sector we consider a class of extended SUGRA 
models with nonuniversalitics in the gaugino masses (see, 
e.g., [HI, [2(3, [2l| and references quoted therein) such 
that the gaugino masses at the scale of grand unification 
are of the form fhi = (1 + 5i)m\/2, i = 1,2,3 for the 
gauge groups U(l) Y , SU(2) L , SU(3) C - Within the above 
framework we discuss illustrative model points: one of 
these leads to an LSP which has a mixed Higgsino-wino 
content while the other is essentially purely wino. It is 
found that the model with the mixed Higgsino-wino con- 
tent can produce the observed positron excess and at the 
same time can coannihilate with the hidden sector gauge 
and chiral fermions to produce a relic density consistent 
with the WMAP [l6| constraints. In contrast, for the 
model where the LSP is essentially pure wino, the above 
phemomcnon fails to bring the relic density within the 
reach of the WMAP data and thus a nonthermal mech- 
anism is needed to generate the right relic density. We 
give now the details of the analysis. 



suppression factors [26 
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Here g A are the degrees of freedom of xa, x = m^o/T x 
with T x the temperature and A A = (m XA — m^o ) / m^o , 
and m^o is the mass of the LSP (x°). The relic abun- 
dance of dark matter at current temperatures obeys the 
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well known proportionality VL^oh 

where xj is the value of x at freeze-out. Because of the 
small couplings of the hidden sector to the visible sec- 
tor, one has <7 aQ , o a p <C a a \,, where a,b = 1, . . . ,n v label 
the MSSM sector particles that are essentially degener- 
ate with the LSP and a, (3 = 1, . . . , rih label the hidden 
sector particles that are also essentially degenerate with 
themselves and the LSP. In this approximation one has 



n^oh 2 



-Bco(^x°/i 2 )mssm, 



(2) 



II. THE GENERAL B Co MECHANISM 

As already noted in the preceding section, a fit to 
the PAMELA data with annihilating dark matter re- 
quires a relatively large annihilation cross section in the 
halo which is as much as two orders of magnitude larger 
(for di-boson production) than the thermal cross section 
needed to fit the data on the relic density of cold dark 
matter (CDM) consistent with WMAP. The main mech- 
anisms discussed in the literature to reconcile the two 
include the enhancement of the velocity averaged annihi- 
lation cross section (av) in the halo cither by annihilation 



near a Breit-Wigner pole 
hancements (see, e.g., |2j 




or by non-pcrturbative en- 
23| and references therein). 

Here we will consider an alternative mechanism con- 
sisting of a set of particles nearly degenerate in mass with 
the LSP but which have suppressed interactions with the 
visible sector states whose particle content is that of the 
minimal supersymmctric extension of the standard model 
(MSSM) . When such a hidden sector is present it can act 
as a reservoir, significantly boosting the resulting thermal 
relic density of the LSP. Here we have in mind a set of 
hidden sector states ac ting as a reservoir, similar in spirit 
to the works of Ref. [Hj] (see also [25[). Additionally 
some states in the MSSM sector may also be degenerate 
in mass with the LSP. In fact for our specific models we 
will find that the light chargino of the MSSM will also 
be degenerate with the LSP. Thus, consider n v number 
of sparticles in the visible sector which arc essentially de- 
generate and rih number of states in the hidden sector 
that are essentially degenerate with the neutralino LSP 
but have suppressed interactions relative to the MSSM 
interactions. In general the relic density is governed by 
°eff = y, ^ b a ABlAlB j where the j A are the Boltzmann 



where (^^o/i 2 )mssm is the relic density as canonically 
calculated using, for example, the standard tools [13] and 
fi^o h 2 is the true relic density when coannihilation effects 
from the hidden sector are taken into account. Further, 
Bco is the enhancement or the boost to the relic density 
that comes from effects of coannihilation in the hidden 
sector. With a aa , o a p <C <J a b , Bq is then given by 
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Here A runs over channels which coannihilate both in the 
MSSM sector and in the hidden sector (i.e., A — l,..,n v + 
rih) and g A are the degrees of freedom for particle A; for 
example, g = 2 for a neutralino and g = 4 for a chargino. 
In the limit that (av) a b are independent of a, b and all 
A A nearly vanish, we find the simple relation 
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Here dh = ^2 a g a is the number of degrees of freedom for 
the coannihilating channels in the hidden sector (with 
suppressed cross sections in the coannihilation process) 
and d v is the number of degrees of freedom in the MSSM 
sector for the coannihilating channels which contribute to 
the coannihilations with the LSP and have interactions 
of normal strength. For the U(l) n hidden sector model 
with each U(l) providing two Majorana states with the 
chargino coannihilating with the LSP, under conditions of 
essentially complete degeneracy of the chargino and the 
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(1 + |n) , while as this degenercy becomes 
the maximal value is Bq q = (1 + 2n) 2 . 



III. SUSY MODELS WITH ENHANCEMENT 
OF THE THERMAL RELIC DENSITY VIA B Co 

As noted above a Bq q of the type discussed in Eq. ((4|) 
can be obtained in a class of models where the hidden 
sector has suppressed interactions with the visible sector. 
Below we construct one explicit example where a hidden 
sector with an extended U(l) n gauge symmetry couples 
with the MSSM sector via the hypercharge gauge field. 
Specifically, we consider an additional contribution to the 
MSSM to be of the form [H, [H 
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where V = {B, X, X', X" . . .} are vector supermulti- 
plets which include the hypercharge gauge multiplet B, 
and $ m are a collection of chiral supermultiplets and 
(Ns,Ny > Ns) are the number of (axions,vectors). 

Thus for the U(V) n extended models we consider Ns — 
n and Ny = n + 1 such that Ns number of vector bosons 
absorb Ns number of axions leading to Ns number of 
massive Z' bosons [H, [29[ . A full analysis would include 
also the electroweak symmetry breaking from the MSSM 
sector generating a mixing between the hidden and the 
visible sectors. Models of this type can arise in a very 
broad class of theories including extensions of the SM, 
extended supergravity models, and in strings and in D- 
brane models. Specifically we will be interested in an 
extension of supergravity models with the hidden sector 
gauge group U(l)x with mixings between the visible and 
the hidden sector given by Eq.((5|). This extension then 
leads to a (4 + 2n) x (4 + 2n) dimensional ncutralino mass 
matrix of the form 
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1 This value of Bq q is the asymptotic limit when there is a com- 
plete degeneracy of the matter in the hidden sector and of the 
chargino with the LSP. There is also an upper asymptotic limit 
which corresponds to a complete split of the chargino from the 
LSP which gives Bq = (l + 2n) 2 , so as the chargino moves from 
a complete degeneracy with the LSP to a complete split, one goes 
from (1 + jj") 2 — ► (1 + 2n) 2 . For n = 3 the above corresponds 
to the transition: Bq = 9 — > Bq q = 49. 



and where a — —cpsy/Mz, P = spswMz, 7 = cpCwMz, 
5 = —spCyyMz- Here (s^c^) = (sin /3, cos (3) , where 
tan/3 = {H-2}/ {H\) and (H2) gives mass to the up quarks 
and (Hi) gives mass to the down quarks and the leptons, 
and (sw,cw) = (sin 6>yi/, cos 6w), where 9w is the weak 
angle. The remaining quantities in Eq.([5]) arc 
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where the ellipses indicate additional U(l)s and where 
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The 4x4 dimensional matrix of Eq.® is the neu- 
tralino mass matrix in the MSSM sector, the (2n) x (2n) 
dimensional mass matrix of Eq. |8]) is the neutralino mass 
matrix in the hidden sector, and the 2n x 4 dimensional 
matrix of Eq. @ is the matrix of off-diagonal terms which 
produce the mixings between the hidden sector and the 
MSSM sector. These mixings arc controlled by the ra- 
tios e = M2/M1, etc. all of which we assume to be much 
smaller than unity [2^|. We also assume that the hidden 
sector neutralinos are all essentially degenerate with the 
LSP neutralino. Thus in the diagonal basis we have the 



set of neutralino states x 



r-0 



1,. . .,4 + 2n) where the 



last rih = 2n are the set of neutralinos from the hidden 
sector. The above serves as an illustrative example, but 
our analysis would apply to any class of models which 
satisfy the conditions discussed in the beginning of this 
section. 

We discuss now two specific model points which fit 
the PAMELA data but lead to significantly different 
sparticle spectra and signatures in the direct detection 
of dark matter and at colliders. 

Higgsino-wino model (HWM) 

Here the soft nonuniversal SUGRA parameters are 

(m ,TOi/ 2 , A)) = (800,558,0) GeV,tan/3 = 5, 

signGu) = +, S h2 . 3 = (-.09, -.50, -.51) , (11) 

and the gaugino-Higgsino content of the LSP is 
X° = 0.726A B - 0.616X W + 0.260/V - O.I6O/12 + C h A x ° h . 
Thus the LSP has a strong Higgsino component in addi- 
tion to strong wino and bino components. Nevertheless 
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we will refer to this model as the Higgsino-wino model 
(HWM) because these components play a major role 
in the analysis to follow. The quantity C/jAx° is the 
component in the hidden sector and is found by excplicit 
calculation to be rather small (\Ch\ < 1% for n = 3). 

Pure wino model (PWM) 

Here the soft parameters arc 

(m ,m 1/2 ,Ao) = (1000,850,0) GeV,tan/3 = 10, 

sign(» = +, $1,2,3 = (0, -.7,0), (12) 

and the gaugino-Higgsino content of the LSP is 
X° = 0.009A B - 0.996Avi/ + 0.081/ii - 0.023/i 2 + C h Ax° h . 
In this model the LSP is dominated by the wino 
component and thus this model will be referred to as 
the (essentially) pure wino model (PWM), while the 
component in the hidden sector C/jAx^ is still small 
(i.e., also \C h \ < 1%, for n = 3). 



there are various cases where boost factors in the range 
~ (1 — 5) fit the PAMELA data in models of supersym- 
metry (SUSY). In Eq. (fT4|) v s is the positron velocity, 
v g ~ c, pq = p(tq) is the local dark matter density in 
the halo (with r Q ~ 8.5 kpc) and lies in the range 
(0.2 - 0.7) GeV/cm 3 [H]. Further, b(E) in Eq.dHJ) is 
given by b(E) = E (E / E ) 2 / t e , with E in GcV and 
Eq = 1 GcV, and where te — t 10 16 s, where r values as 
large (1-5) have been considered in the literature and we 
adopt t = 3 [33| . Here (av) h a io is the velocity averaged 
cross section in the halo of the galaxy and X(E,E') is 
the halo function. We have considered both the Navarro, 
Frenk and White (NFW) and Moore et. al [HJ profiles 
coupled with various diffusion models. The diffusion loss 
equation depends on the propagation parameters 8, Kq 
solved in the region modeled with cylindrical symmetry 
bounding the galactic plane with height 2L. The dimcn- 
sionless halo function has been parameterized to satisfy 
constraints on the boron to carbon ratio 0, [3l| such that 



For both model points given above one finds the fol- 
lowing mass relations 



A A2 



X-3 



A4 



(13) 



These mass relations follow simply from the condition 
that rh 2 < fh\, \fi\ >> Mz- However, a distinguishing 
feature of the HWM is that the mass gap between the 
LSP and the chargino is order ~ 10 GeV while for the 
case of the PWM the mass gap is order the pion mass. 
The neutralino annihilation in both models is dominated 
by W + W~ production in the halo. We will discuss this 
in much more detail in what follows. 
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and where £ = log 10 A d with A^i in units of kpc 
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IV. THE PAMELA POSITRON DATA AND 
THERMAL NEUTRALINO DARK MATTER 

One issue of interest in this work is to address the 
compatibility of the recent PAMELA positron data and 
the WMAP data. Thus, we begin with a brief review of 
the calculation of the primary positron flux relevant for 
the models considered here. 

The positron flux which enters as a solution to the 
diffusion loss equation under steady state conditions is 
given by [H,^,! 



ME) = 



8nb(E) m\ 



^-F(E) , 



(14) 



F(E) = 



r M-o 
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In the above, B s is a so-called boost factor which pa- 
rameterizes the possible local inhomogeneities of the 
dark matter distribution. Large boost factors have been 
used in the literature, even as large as 10,000, to ex- 
plain the PAMELA data. However, we will show that 
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MIN(M2) 0.55 



0.00595 
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wf 




wl 


- 2.28838 


-0.605364 


- 0.287614 


-0.762714 


w% 


tug 






- 0.319561 


-0.0583274 


- 0.00503555 


-0.00016691 



TABLE I: The class of models consistent with the 
boron/carbon ratio [3lJ considered here. The MED and MAX 
halo models are not shown since they tend to overproduce 
the positron flux at low energies. Also shown are the co- 
efficients (a,i,bi) entering in the dimensionless halo function 
for positron propagation (j§] where the NFW profile is used 
and the MIN(M2) diffusion parameters enter. The last set of 
entries are the fragmentation coefficients wf for the W + W~ 
mode 13011. 



For the case of dark matter in the halo which primarily 
annihilates into W boson pairs, the positron fragmenta- 
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FIG. 1: SUSY predictions for the positron flux ratio: The top 
curve is for the PWM (see Eq. (JT3J) where the neutralino is 
wino dominated, with m^o = 199 GeV, and (av)ww — 1.96 x 
10~ 24 cm 3 /s with a Z? c iump of 1.0 and Q.h 2 ~ .01 with a Bco ~ 
9 (the asymptotic value for the wino case). The middle curve 
is for the HWM (see Eq. (jXTJ ) where the neutralino has a 
mixed wino, bino and Higgsino content, with m^o = 195 GeV, 
and {av)ww = 0.28 x 10 _24 cm 3 /s with B c ium P of 4 and gives 
rise to an fi/i 2 = .085, with a Bco of 16. We have taken 
Pq = 0.6 GeV/cm 3 and r = 3 for both curves. Slightly larger 
clump factors can easily accommodate a downward shift in 
the product P%tb- The bottom curve is the background. The 
experimental data is from [l[. 



tion function can be expressed as [3C 



dm 
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exp 
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where x = E/m and the fragmentation function was fit 
in the analysis of [H with HERWIG [H. 

The positron fluxes in the absence of dark matter anni- 
hilations have been parameterized in [36j using fits to the 
analysis which we adopt here. For the case of positrons 
annihilating strongly into W bosons, we find that the 
MED and MAX models with NFW or Moore profile tend 
to overproduce the positron flux data at energies well 
above the region of solar modulation. Thus the class of 
models we consider restricts the profile/diffusion model 
to a MIN(M2) scenario, at least for the case of positrons. 
We thus display only the fit parameters for this scenario 
in Table jT|. The fit parameters for the halo function 
are also shown in Table jlj along with the fragmenta- 
tion functions for dark matter annihilations into W + W~ 
states. As discussed in the previous section, in the anal- 
ysis here we fix the soft parameters of the SUGRA models 
at the scale of grand unification and examine the implica- 
tions of the model at the electroweak scale. The analysis 
shows that the mixed Higgsino-wino model (HWM) can 



FIG. 2: An illustration of the Bco-mechanism with a con- 
tour plot of constant Bco in the Ayisibic vs AHidden plane for 
the U(l) 3 extended SUGRA model. Large values of Bco are 
admissible for the mixed Higgsino-wino case (see the text). 
The region to the right of the contour Bco = 9 is the domain 
of the PWM while the region to the left of it is the domain 
of the HWM. Thus the largest value of B Co in the PWM is 9 
while for the HWM case Bco can be significantly larger. 



have mass splittings between the LSP and the chargino 
of around 10 GeV and still produce relatively large halo 
cross sections that do not require large boost factors to fit 
the PAMELA data. On the other hand for the pure (or 
nearly pure) wino LSP one finds that the mass splitting 
of the LSP neutralino and the light chargino is negligi- 
bly small and this result is generally a necessary conse- 
quence of the wino dominated class of models. A numer- 
ical analysis of the positron excess for comparison with 
the PAMELA experiment is presented in Fig. ([T]) for the 
HWM and for the PWM. It is found that both models 
are capable of fitting the PAMELA data. A wino LSP of 
a mass near 200 GeV has been suggested in Ref. [|[ as a 
good candidate for explaining the positron excess using 
a methodology and choice of parameters similar to the 
ones used here, though the model in that study and this 
study are quite different. 

While both the HWM and PWM are successful at 
reproducing the PAMELA positron excess, the thermal 
relic densities for the two models are very different. For 
the case of the HWM, the relic density can be consistent 
with the WMAP data for a thermal LSP. However, the 
relic density for the PWM from thermal processes is at 
best about a factor of 10 too small and here one needs a 
non-thermal mechanism to be compatible with WMAP, 
as was assumed in the work of Ref. [|[ . The compatibil- 
ity of the relic density with WMAP for the HWM comes 
about because of the enhancement from the hidden sec- 
tor, i.e., because of the Bc factor. As seen in Fig. |(3J), 
Bqo can acquire values as large as 45. This limiting value 
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Model 8 K (kpc 2 /Myr) L (kpc) V mm (km/s) 



10 



MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.70 


0.0112 


4 


12 


MAX 


0.46 


0.0765 


15 
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TABLE II: Specific models considered here, consistent with 
the boron/carbon ratio [3f§] 



of Bco is realized only when there is a complete split be- 
tween the LSP and the light chargino. While this value 
is not achieved in the HWM, the presence of a ~ 10 GeV 
split between the LSP and the chargino does conspire to 
produce a Bq q in the range 16 — 20 while for the PWM 
the maximal value one may have is Bq ~ 9. While a fac- 
tor of 16 — 20 is good enough to boost the relic density for 
the HWM to be consistent with WMAP, for the PWM 
the relic density still falls below the WMAP corridor and 
would require a very large collection of degenerate states 
(e.g., n = 11) to bring the relic density within the corri- 
dor. 



V. ANALYSIS OF THE p FLUX 

Quite similar to the analysis of the positron flux, the 
anti-proton flux assumes a steady state solution to a dif- 
fusion equation, retaining a factorizcd form split between 
the particle physics content and the astrophysical details, 
the latter encoding the halo and the propagation model. 
Thus the anti-proton flux from annihilating ncutralinos 
in the galaxy is given by [I?) 



%(T) = B, 



P St: m| 



R(T) Y,' 
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halo 



~dr~ 



(20) 



where T is the kinetic energy, and R(T) has been fit as in 
Ref. 0] for various profile/diffusion models. Specifically 
for the NFW models one has 

MIN : log 10 [i?(T)/Myr] = 0.913 + 0.601 r 
-0.309 r 2 - 0.036 r 3 + 0.0122 r 4 

MED : log 10 [i?(T)/Myr] = 1.860 + 0.517 r 
-0.293 r 2 - 0.0089 r 3 + 0.0070 r 4 

MAX : log 10 [i?(T)/Myr] = 2.740 - 0.127 r 

-0.113r 2 + 0.0169 r 3 - 0.0009 r 4 , (21) 

where r = r(T) = log 10 (T/GeV), and where the p prop- 
agation parameters are given in Table dTTJ) . 

As discussed already, for the models of interest we 
study here, the contribution to (av) in the halo is dom- 
inated by several orders of magnitude from the annihi- 
lations of ncutralino dark matter into W + W~ bosons. 
The p fragmentation functions needed for the analysis of 
the flux are given by [39j|,[3(| as 



dN, 



ww 



dx 



(Pia^+pal logical)- 1 , (22) 
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FIG. 3: The absolute anti-proton flux for the HWM with var- 
ious halo/diffusion models NFW(MIN,MED,MAX) and only 
(MIN, MED) for the PWM . Here p = 0.3 GeV/cm 3 (the sig- 
nal scales as p 2 ), Bp = 1. For the background flux (band) we 
have adopted the parameterizations of Ref. [4(j. Also shown 
is the preliminary PAMELA data [n], along with the ear- 
lier data sets from BESS and CAPRICE The analysis 
here shows that the models discussed in this work can accom- 
modate the anti-proton flux constraints. 



where here x = T/m^o and T is the kinetic energy The 
parameters pi in the above equation depend on the neu- 
tralino mass and are given by 



Pi (m) = (a. a m a ' 2 + a a m a ^) 1 , 
where the values of are given in Table (III) 



m = m^o , (23) 





3 = 1 


3=2 


J = 3 


J = 4 


i — 


1 


306.0 


0.28 


7.2xl0" 4 


2.25 


i — 


2 


2.32 


0.05 








i — 


3 


-8.5 


-0.31 








i — 


4 


-0.39 


-0.17 


-2.0 x 10" 2 


0.23 



TABLE III: Coefficients CHj for W + W pro cess entering into 
the anti-proton fragmentation functions [39l| [30l] . 



The anti-proton flux observed at the top of the atmo- 
sphere including solar modulation can be accounted for 
by replacing $ P (T) = $p(T,0 F ) = $ p - (T + \Z\<f> F ) and 
including a kinetic energy correction ratio, such that 



$®(T) = 



((T + m p ) 2 - ml) ■ $g (T + \Z\</> F ) 
(T+\Z\4> F + m p y-ml 



(24) 



where m p is the proton mass, Z = 1, and the Fisk po- 
tential <j) F is taken as 500 MV. In Fig.Q we give an 
analysis of the anti-proton flux with the parameters in- 
dicated in the caption of the figure showing both the 
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SUSY signal and the backgrounds. The analysis demon- 
strates that HWM is only currently constrained for the 
MAX diffusion model, while the PWM is constrained for 
a MED diffusion model but is essentially unconstrained 
for a M1N diffusion model. Similar conclusions regarding 
the sensitivity of the p flux to the halo/diffusion model 
have also been made in the first Ref. of Q. 

In summary, the findings here for the PWM are gen- 
erally in good agreement with the recent results of [J], 
and further we find that the HWM is only very weakly 
constrained by the PAMELA anti-proton data as well as 
by the earlier BESS and CAPRICE data. 



VI. PHOTON FLUX; EGRET AND FERMI-LAT 




E GeV 



The continuum photon source from the annihilating 
SUSY dark matter is given by [iij 



d 2 $ 7 
dn dE 



1 r_Q p% 

2 47rm| 



E 



(TV, 



halo 



dN* 

-di 



(25) 



l 2 +r 2 , 



where J{ip) = (r Q p|) 1 J °° p 2 (r)dl(ip) and r 2 - t -p, - 
2Zr0Cos?/> with ip the angle of integration over the line 
of sight. After integration, the astrophysics is encoded 
in J • A£7, where J = (1/Afi) J An </(?/>) d£l and various 
values of J are given for a collection of angular maps, 
in, for example [45[ (see also [III, [I?) , [H[ for early work) . 
As remarked previously, W boson pair production is the 
dominant contributor to the photon flux in the models wc 
discuss, and an upgraded set of fragmentation functions 
are given in Q and are shown below 



E- 



dE 



exp 



\ - w 



^-hx n (E/M) 



(26) 



Separately, effects due to Brehmstrahlung have been 
studied in Ref. [49[ and arc easily accounted for. Val- 
ues of relevant for the analysis are listed below. 



-6.751 -5.741 -3.514 -1.964 -0.8783 -0.2512 -0.03369 



FIG. 4: An exhibition of the gamma ray flux the HWM and 
the PWM in the coordinate range indicated, with the Einasto, 
NFW, and Isothermal profiles with the same halo cross sec- 
tions as given in Figs. (|1I3 |) . Shown are the EGRET Q results 
and FERMI-LAT results as reported in [jy], H2 along with 
the background flux (band). The analysis here shows that 
both models can accommodate the photon flux constraint. 



between the HWM and the PWM. Here the PWM will 
produce a significantly stronger monochromatic source 
than the HWM which is dictated by the size of the rela- 
tive cross sections since m^o are essentially the same for 
the HWM and PWM: 

HWM : M^ loop = 1.6 x 1CT 28 cm 3 /s, (27) 



PWM 



^-^p = 2.0 x 10~ 27 cm 3 /s. 



™)^ l0 ° P = 1-0 x 10_27 cm 7 s < (28) 



PWM : (crv)l /Z loop = 1.3 x KT 26 cm 3 /s. 

Thus while the photon energies will be essentially the 
same (m^o for the 77 channel and m^o(l — S) , S = 
M|/(2m^o) 2 for the 7Z channel) in fact, the PWM 
would predict a flux an order of magnitude larger than 
the HWM in both channels. 



Our result for the photon flux is given in Figure ([4]) where 
we exhibit the gamma ray flux arising from the annihi- 
lation of neutralino dark matter in the 10 — 20 region 
(lOdeg < \b\ < 20deg), (Odeg < I < 360 deg) [latitude b 
and longitude I] which is the region relevant for compar- 
ison with the preliminary FERMI-LAT results [EH . Fig- 
ure ([3]) shows that the more accurate FERMI-LAT data 
falls in magnitude below the EGRET data. Our anal- 
ysis of the photon flux for the HWM is consistent with 
the FERMI-LAT and EGRET data, while the PWM may 
show an excess at extended energy ranges. Note the max- 
imum of the flux for the signal appears at the last data 
point of the FERMI data near 10 GeV. 

We note that monochromatic sources [HH (calculated 
with DarkSusy [27| ) yield a further distinguishing feature 



VII. EFFECTS ON DIRECT DETECTION 

The direct detection of dark matter is sensitive to the 
Higgsino content of the LSP. For the HWM the LSP has 
a significant Higgsino component and thus the spin in- 
dependent cross section for neutralino-proton scattering 
in the direct detection dark matter experiments is sig- 
nificant and may lie within reach of the next generation 
experiments. For the PWM, we remind the reader that 
the LSP is essentially 100% wino, and in this case the spin 
independent cross section will be very small, essentially 
beyond the limit of sensitivity of near future experiments 
on the direct detection of dark matter. Specifically, a di- 
rect calculation using MicrOmegas (see Ref. 2 of [13|) 
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yields a ncutralino-proton cross section of 

HWM : crp p = 6.62 x l(T 8 pb, (29) 
PWM : af 0p = 1.15 x 10~ 9 pb, (30) 

which translates into a recoil rate on germanium targets 
of 

HWM: R= 1.58 X 10" 2 /day/kg, (31) 
PWM: R = 2.81 x l(T 4 /day/kg, (32) 

and on liquid xenon targets of 

HWM: R = 2.60 x l(T 2 /clay/kg , (33) 
PWM : R = 4.63 x 10~ 4 /day/kg j (34) 

when integrated over the entire range of nuclear recoil 
energies. For the pure wino and pure higgsino cases the 
cross sections can receive important loop corrections [55]. 
From the analysis of [55| one can estimate these effects, 
and they are found not be significant for the parameter 
space we investigate. Clearly, the prospects for the di- 
rect detection of the relic LSP are most promising for 
the HWM, which should be accessible in near future 
dark matter experiments with an improvement in sen- 
sitivity by a small factor, as is expected to occur. For 
the PWM the observation of the spin independent cross 
section requires an improvement in sensitivity by a factor 
of about 100 which is unlikely to happen in the foresee- 
able future. 



VIII. COLLIDER IMPLICATIONS OF PAMELA 
INSPIRED MODELS 

Discovery Modes at the LHC: An important and ro- 
bust aspect of supersymmetric models which are capable 
of generating the observed PAMELA positron excess is 
that some supcrpartners will necessarily be light and are 
therefore potentially discoverable at the LHC. Certain 
key physical masses for the HWM and PWM are given in 
Table pvp . All masses are computed from the high-scale 



Mass 


HWM 


PWM 




Mass 


HWM 


PWM 




198.9 


195.2 






648.5 


1516 




217.0 


357.0 




m i 2 


866.8 


1749 


x 3 


429.9 


1025 






841.4 


1729 


th~o 


451.3 


1029 






970.2 


1902 


m~± 


208.8 


195.5 




rrif 1 


817.7 


1011 


m~± 
m-g 


448.6 


1036 




m? 2 


822.8 


1041 


707.1 


1929 











TABLE IV: Relevant SUSY mass spectra for the HWM and 
PWM as calculated from the high-scale boundary conditions 
given in flllj l and (|12p . respectively. All masses are in GeV. 



boundary conditions of Eq. (fTTj) and Eq. (fl"2|) after renor- 
malization group evolution using SoftSUSY (54[. One 
finds that for the PWM the sum rules of Eq. (fl~3|) are sat- 
isfied with an accuracy of less than 0.5%. In fact, for 
the PWM there is an almost perfect degeneracy between 
the lightest chargino and the lightest neutralino mass, as 
is typical for models with a wino-dominated LSP. In the 
HWM this mass difference is larger, reflecting the larger 
proportion of bino and Higgsino components for the LSP 
wave function, so that Eq. (fT3")) is satisfied only at the 5- 
6% level. But as we saw in Section [TV] this results in a 
larger Bqo factor for the HWM. 

One may note that in Table pV[) a significant differ- 
ence in mass scales for the SU (3)-charged supcrpartners 
of the HWM and of the PWM. This difference in the mass 
scales has large implications for the discovery prospects 
of the two models. To analyze the signatures of these 
models at the LHC we generated events using PYTHIA 
followed by a detector simulation using PGS4 [56|. Two 
data sets for each model at ^fs = 14 TcV were gener- 
ated for 10 fb _1 and 100 fb" 1 of signal events, as well 
as a 500 pb _1 sample at y^s = 10 TeV for each model. 
In addition we considered a suitably-weighted sample of 
5 fb _1 Standard Model background events, consisting 
of Drell-Yan, QCD dijet, ti, bb, W/Z+jets and dibo- 
son production. Events were analyzed using level one 
(LI) triggers in PGS4, designed to mimic the CMS trig- 
ger tables d^j. Object-level post-trigger cuts were also 
imposed. We require all photons, electrons, muons and 
taus to have transverse momentum pt > 10 GeV and 
\rj\ < 2.4 and we require hadronic jets to satisfy \rj\ < 3. 
Additional post-trigger level cuts were implemented for 
specific analyses, as described below. 

We begin with standard SUSY discovery modes [Hj], 
slightly modified to maximize the signal significance for 
these models. These five signatures are collected in Ta- 
ble (|V| for 10 fb _1 of integrated luminosity. In all cases 
we require transverse sphericity St > 0.2 and at least 
250 GeV of $t except for the trilepton signature, where 
we place a cut of Pt > 200 GeV. The multijet chan- 
nel includes a veto on isolated leptons and requires at 



HWM PWM 



Signature 


Events 


s/Vb 


Events 


s/Vb 


Multijets 


8766 


183.74 


50 


1.05 


Lepton + jets 


2450 


32.25 


26 


0.34 


OS dileptons + jets 


110 


6.39 


4 


0.23 


SS dileptons + jets 


60 


11.77 





NA 


Trileptons + jets 


14 


2.47 





NA 



TABLE V: LHC discovery channels for the HWM and the 
PWM: Event counts are after 10 fb _1 of integrated lumi- 
nosity. All signatures require transverse sphericity St > 
0.2 and at least 250 GeV of Pt except for the trilepton 
signature, where only JEt > 200 GeV is required. Here 
(OS,SS) = (opposite sign, same sign) 
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IVL, Distribution 




Effective Mass (GeV) 

FIG. 5: Effective mass distribution for the Higgsino-wino 
mixed model(HWM) for 5 fb^ 1 (yellow) and the wino model 
for 100 fb _1 (blue) along with the Standard Model back- 
ground (dashed open histogram). 



least four jets with the transverse momenta of the four 
leading jets satisfying pr > (200, 150, 50, 50) GeV, re- 
spectively. For the lcptonic signatures we include only 
and /i^ final states and demand at least two jets with 
the leading jets satisfying pr > (100,50) GeV, respec- 
tively. For the case of the HWM, the existence of light 
squarks and gluinos give rise to a relatively large num- 
ber of events passing the cuts, and thus the model is 
easily accessible at the LHC in almost all channels. In 
fact, wc find that this particular model point gives rise 
to ~ 150 multijet events with $t > 250 GeV with just 
500 pb^ 1 at yfs = 10 TeV. By contrast the PWM has 
nearly 2 TeV squarks and gluinos and requires over a 100 
fold increase in luminoisty to reach a comparable event 
rate in this channel. Thus after 100 fb _1 of integrated 
luminosity the PWM results in only 440 multijet events 
with E T > 250 GeV . This is illustrated in Figure © 
where the effective mass, defined as the scalar sum of the 
transverse momenta of the four hardest jets in the event 
plus the missing transverse energy, is plotted for events 
satisfying the cuts described above for the multijet chan- 
nel. The signal for the HWM is clearly discernible above 
the Standard Model background after 5 fb _1 . For com- 
parison we plot the same distribution for the PWM after 
100 fb _1 . 

Experimental Challenges of the PWM Model: The sup- 
pression of leptonic final states for both models, espe- 
cially the trilepton channel, is the result of the small 
mass difference between the low-lying electroweak gaug- 
ino states which causes leptonic decay products to be gen- 
erally quite soft. This is particularly severe for the wino- 
like scenario (the PWM) . We note that the total (leading 
order) supersymmetric cross section for the PWM is still 



HWM PWM 



Object Cuts (GeV) 


Events 


s/Vb 


Events 


s/Vb 


p^ > 150, #t > 150 
P>t > 200, $ T > 150 


1994 
1302 


2.13 
2.52 


3442 
1983 


3.68 
3.84 


p£* > 150, # T > 200 
p£* > 200, # T > 200 


1334 
1241 


2.53 
2.58 


2147 
1904 


4.08 
3.95 


p'r > 150, # T > 300 


659 


3.57 


771 


4.17 



TABLE VI: Monojet signature for the HWM and the PWM: 
Event counts are after 100 fb _1 of integrated luminosity for 
various choices of cuts on total $t and jet pr- All signatures 
involve a lepton veto and require no other jets in the event 
with pJp > 20 GeV. No transverse sphericity cut was applied 
in any of these signatures. 



a healthy 2.3 pb (to be compared with 7.4 pb for the 
HWM). The lack of a signal here is largely the result 
of an inability to trigger on events in which light elec- 
troweak gauginos are produced (99% of the total SUSY 
production cross-section) and the absence of hard leptons 
in the decay products of these states. These difficulties 
are common to phenomcnological studies of models where 
the mass gaps between sparticles are small [2(| [59|, HO, Hl| 
such as in models with anomaly-mediated supersymmc- 
try breaking (AMSB) [62j which share many of the same 
features in the gaugino sector [g3, H3 as the PWM. 

Nevertheless, some signatures unique to this scenario 
can be explored with sufficient integrated luminosity. For 
example, one can look for events which pass the initial 1£t 
trigger, but which have no leptons or energetic jets. Such 
events predominantly arise from direct production of Xi 
and/or Xi 2 states whose decay products produce only 
very soft objects. If the mass difference Am = m~± — m^a 
is sufficiently small the chargino will travel a macroscopic 
distance before decaying. An offline analysis may then re- 
veal events with tracks in the inner layers of the detector 
which abruptly end, lea ving no further leptonic tracks or 
calorimeter activity [t35l . l66l . [6^ . In the extreme limit of 
a pure-wino LSP, as in the AMSB models, this mass dif- 
ference is typically the size of the pion mass, i.e. OijnV) 
and the distance traveled can be several centimeters. In 
the case of the PWM, however, the mass difference is 
roughly twice the pion mass and the typical decay length 
will be such that the decays will appear prompt and few 
events will reveal a displaced vertex. More promising 
is to consider events with a single high-p^ jet and large 
missing transverse energy. Such events can arise from 
initial state radiation in electroweak gaugino production, 
or in cases where the lightest chargino or neutralino is 
produced in association with a gluino or squark. This 
particular 'monojet' channel has relatively large event 
rates for both the HWM and PWM, with the primary 
Standard Model background coming from TV+jets pro- 
duction. Event rates and signal significance for various 
jet px and $t cuts are given in Table (|VI[) . More detailed 
analyses of similar models show that such events can be 
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Opposite-Sign Dilepton Invariant Mass 
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Two b Jets Invariant Mass 



600 800 



1000 



5120 



Model HWM: 100 fb 



Same Flavor 



Flavor-Subtracted 




FIG. 6: The invariant mass distribution of 2 b-jets events 
for HWM where the LSP contains substantial Higgsino and 
wino components. Shaded histograms represents the Stan- 
dard Model backgrounds. The solid lines are the best fit 
functions for the SUSY signal near the endpoint. The end- 
point is estimated to be ~ 530 GeV based on the linear fitting 
functions. The embedded window shows the mass distribu- 
tions for both SUSY and SM, when one requires a 200 GeV 
missing energy cut. In order to suppress the Standard Model 
background, we increase the missing energy cut to 400 GeV 
and add two additional conditions: (1) lepton veto; (2) at 
least 2 more jets besides the 2 b-tagged jets. 



bona-fidc discovery modes in scenarios of this type [681 ] . 

Measuring Sparticle Masses in the HWM Model: For 
the mixed Higgsino- wino model (HWM) a sufficient num- 
ber of events can be obtained even with a low integrated 
luminosity which will allow one to determine the prop- 
erties of the superpartner spectrum and can confirm the 
features of Table pV[) - particularly those with direct 
relevance to the calculation of the thermal relic abun- 
dance and positron yield in cosmic rays. Thus even 
with 10 fb -1 it should be possible to get a reasonable 
estimate of the gluino mass by considering events with 
precisely two b-tagged jets. As the gluino is relatively 
light it would be produced in significant amounts and 
the invariant mass distribution of b-jets produced in its 
three-body decays will reveal a kink which allows one 
to determine the lower limit on the gluino mass know- 
ing the LSP mass (69j . i.e. m„ > 



, kink 



Over 3000 events with two b-jets satisfying p 1 ^ > 60 
GeV and $ T > 200 GeV were produced in our 10 fb" 1 
sample for the HWM. With these cuts, the Standard 
Model background, arising mostly from ti events, is com- 
parable to the signal as is clear from the inset of Fig- 
ure ([6]). To reduce this background we veto events with 
isolated lcptons, require two additional jets without b- 
tags, each satisfying p!^, > 60, and increase the missing 
transverse energy cut from 200 GeV to 400 GeV. This 



FIG. 7: Dilepton invariant mass distributions for the HWM 
after 100 fb -1 . The unshaded histogram gives the invariant 
mass distribution for events with precisely two opposite-sign 
leptons of the same flavor (electron or muon). The shaded 
histogram gives the flavor-subtracted distribution which re- 
sults when the opposite-sign, opposite-flavor distribution is 
subtracted from the same-flavor distribution. 



reduces the signal sample by approximately a factor of 
four, but the Standard Model background is now reduced 
to manageable levels. This is displayed in the main panel 

of Figure © where (Mg v ) kink ~ 500 GeV, which gives 
m g <; 700 GeV, which is consistent with the true gluino 
mass of 707 GeV in this model. 



Of even more relevance is the relatively small mass gap 
between the two lightest neutralinos \2 an d X ■ While 
not as severe as in the wino-like case of the PWM, the 
mass difference is still small enough that the number of 
energetic leptons coming from the final stages of spar- 
ticle cascade decays will be suppressed relative to the 
number expected from more universal models. The re- 
duction in events with two or more energetic, isolated 
leptons will significantly degrade the ability to make mea- 
surements of mass differences using the edges of various 
kinematic distributions. For example, a typical strat- 
egy for accessing the mass difference between light neu- 
tralinos is to form the flavor-subtracted dilepton invari- 
ant mass for events with at least two jets satisfying 
Pj? > 60 GeV, at least 200 GeV of $t and two opposite- 
sign leptons [70| | . The invariant mass distribution formed 
from the subset involving two leptons of opposite flavor 
is subtracted from that involving two of the same flavor, 
i.e., the combination (e + e~ + — e + fj,~ — e~/i + ), 

to reduce background. In Fig. J7J) we plot the invariant 
mass of same-flavor, opposite-sign leptons in two-lepton 
events for HWM after 100 fb -1 as well as the flavor- 
subtracted distribution. Note the small number of events 
which remain after the subtraction procedure has been 
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performed. Nevertheless, the beginnings of a feature in 
the low-energy bins can be discerned which is consistent 
with the mass differences between the low-lying gaugi- 
nos in this model. Additional statistics and more careful 
analysis techniques will be needed to make strong state- 
ments about the neutralino masses in this model. 

ILC Implications: The ideal machine to study the spec- 
troscopy of the light gauginos and to confirm the model 
as a potential explanation for the PAMELA positron ex- 
cess would be at a potential International Linear Collider 
(ILC). We emphasize that a linear collider operating at 
a center-of-mass energy of 500 GeV would be sufficient 
to study the closely-spaced lightest neutralinos and light- 
est charginos of either model. For the case of the pure 
wino model (PWM) a future linear collider will be essen- 
tial for resolving the presence of two nearly-degenerate 
states near 200 GeV and for studying their couplings. 
The prospects for both models at the ILC are quite good 
with a{e + e- -> SUSY) = O(0. 1) pb for both models. Fi- 
nally, it is worth pointing out that for both models there 
is the distinct possibility of observing a degenerate clus- 
ter of Z' bosons with masses of the order of the lightest 
chargino mass. This will only be possible if the decays 
of such additional Z' bosons into hidden sector matter 
states are forbidden or largely suppressed. In such cases 
they would appear as sharp resonances in Drell-Yan pro- 
cesses and should be discernible at the LHC with suffi- 
cient luminosity. Observation of such states would be a 
spectacular confirmation of the Bq mechanism of gen- 
erating the correct neutralino relic density by thermal 
means for the HWM model. 

IX. CONCLUSIONS 

The positron excess observed in the PAMELA satel- 
lite experiment has spawned various mechanisms to ex- 
plain this effect. An interesting possibility relates to 
the positron excess arising from the annihilation of dark 
matter in the galaxy. An adequate explanation of this 
phenomenon within a particle physics model then would 
require a simultaneous fit both to the WMAP data re- 
garding the density of cold dark matter as well as to 
the PAMELA data. Often it turns out that models that 
give the desired relic density give too small a (av) in 
the galaxy to explain the PAMELA excess. Alternately, 
models which produce an adequate (av) and explain the 



PAMELA excess fail to produce the proper relic den- 
sity. To reconcile the two phenomena typically the fol- 
lowing mechanisms have been proposed: (i) For models 
that generate the right relic density but give too small a 
(av) the Sommerfeld enhancement or the Breit-Wigner 
pole enhancement of (av) can explain both sets of data; 
(ii) For models which give an adequate (av) to explain 
the PAMELA excess but give too small a relic density, a 
non-thermal mechanism is taken to produce the proper 
relic abundance. 

In this paper we have carried out an analysis of a 
distinctly different mechanism which can generate the 
proper relic density for the second type of models, 
i.e., models where (av) is large enough to explain the 
PAMELA excess but the relic density needs an enhance- 
ment. We illustrate this mechanism within the frame- 
work of nonuniversal SUGRA models with an extended 
hidden U(l)\ gauge symmetry. The extension along 
with conditions of mass degeneracy in the hidden sec- 
tor gives rise to predictions which match the PAMELA 
data and the predicted relic density is consistent with the 
WMAP data with a mixed Higgsino-wino LSP. The anti- 
proton and the gamma ray fluxes emanating from the 
annihilation of dark matter in the galaxy are found to be 
compatible with data. Implications of the model for the 
direct detection of dark matter and some of its collider 
signatures were also discussed. The model is testable on 
both fronts. Specifically the mixed Higgsino-wino LSP 
model is testable at the LHC with just 1 fb _1 of data. 
While the analysis was done in the framework of nonuni- 
versal SUGRA models, the results of the analysis are 
valid in a broad class of models including string and D- 
brane models as long there is an LSP with mass in the 
proper range and a suitable admixture of Higgsino-wino 
components with an appropriate degeneracy in the hid- 
den sector. Finally, we note that we have not made an 
attempt here to fit the FERMI-LAT e+ + e~ result 0. 
Such a fit can be accommodated by assuming that the 
high energy flux is a consequence of pulsars or mixed 
dark matter and pulsar contributions [H, [7l| . 
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